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Abstract: Naturalness demands that the quadratic divergence of the one-loop top 
contribution to the Higgs mass be cancelled at a scale below 1 TeV. This can be achieved 
by introducing a fermionic (spin- 1/2) top partner, as in, for example, Little Higgs 
models. In this paper, we study the phenomenology of a simple model realizing this 
mechanism. We present the current bounds on the model from precision electroweak 
fits, flavor physics, and direct searches at the LHC. The lower bound on the top partner 
mass from precision electroweak data is approximately 500 GeV, while the LHC bound 
with 5 fb _1 of data at y^i = 7 TeV is about 450 GeV. Given these bounds, the model 
can incorporate a 125 GeV Higgs with minimal fine-tuning of about 20%. We conclude 
that natural electroweak symmetry breaking with a fermionic top partner remains a 
viable possibility. We also compute the Higgs decay rates into gauge bosons, and find 
that significant, potentially observable deviations from the Standard Model predictions 
may occur. 
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1 Introduction 

The Standard Model (SM) of particle physics postulates the existence of an elementary 
scalar field, the Higgs, which is responsible for electroweak symmetry breaking. Pre- 
cision measurements of the properties of electroweak gauge bosons are consistent with 
this picture, and favor a light (~ 100 GeV) Higgs boson. Recently, experiments at the 
Large Hadron Collider (LHC) reported preliminary evidence for a new particle with 
properties roughly consistent with the SM Higgs and a mass of about 125 GeV [1]. 

In the SM, the contribution of quantum loops to the Higgs mass term is quadrat- 
ically divergent. To avoid fine-tuning, new physics beyond the SM must appear and 
cut off this divergence at a scale of order 1 TeV or below. Precision electroweak data 
favors models where the divergence is cancelled by loops of new weakly-coupled states; 
such cancellations can occur naturally as a consequence of underlying symmetries of 
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Figure 1. One- loop Higgs mass renormalization in a model with a fermionic top partner, 
such as the Littlest Higgs. 



the theory. What is the minimal set of new particles that must appear below 1 TeV to 
avoid fine-tuning? It is well known that the only SM contribution to the Higgs mass 
that must be modified at sub- TeV scales is the one-loop correction from the top sector. 
All other SM loops are numerically suppressed by either gauge or non-top Yukawa cou- 
plings, by extra loop factors, or both. As a result, the states responsible for cutting off 
these loops can lie above 1 TeV with no loss of naturalness. Thus, the sub- TeV particles 
that soften the divergence in the top loop, the "top partners," provide a uniquely well- 
motivated target for searches at the LHC, and it must be ensured that a comprehensive, 
careful search for such partners is conducted. 

The best-known mechanism for canceling the Higgs mass divergences is super- 
symmetry (SUSY). In SUSY models, the quadratic divergence in the SM top loop is 
cancelled by loops of scalar tops, or stops. Recently, a number of papers [2] empha- 
sized the importance of stop searches at the LHC, and reinterpreted the published 
LHC results, based on the 1 fb _1 integrated luminosity data set, in terms of bounds 
on stop masses. It was found that completely natural spectra are allowed so far. On 
the other hand, incorporating a 125 GeV Higgs in the Minimal SUSY Model (MSSM) 
does require significant fine-tuning, of order 1% at best. (Fine-tuning can be reduced 
in non-minimal models [3].) 

However, SUSY is not the only option for canceling the quadratic divergence in the 
SM top loop. An alternative is to introduce a spin-1/2 top partner T, a Dirac fermion 
with mass ro^, which is an SU(2)l singlet, color triplet, and has electric charge 2/3. 
In the Weyl basis, T = (Xl,Tr). This field couples to the SM Higgs doublet H via 

£ = -XtT^HQs + %^{H^H)TlT R + h.c. , (1.1) 

where Q3 is the SM third-generation left-handed quark doublet, Xt is the SM top 
Yukawa, Ay is a new dimensionless coupling constant, and H = {ia2Hy. The one- 
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loop contribution to the Higgs mass in this model is shown in Fig. 1; the quadratic 
divergences present in each of the three diagrams cancel in the sum. Even though the 
structure of the couplings in Eq. (1.1) looks completely ad hoc at first sight, it can 
emerge naturally if the Higgs is embedded as a pseudo-Nambu-Goldstone boson [4] of 
spontaneous global symmetry breaking at the TeV scale. The global symmetry must 
be broken explicitly to induce non-derivative Yukawa and gauge couplings of the Higgs; 
divergence cancellation is achieved if the explicit symmetry breaking terms obey the 
"collective" condition, such as in Little Higgs models [5, 6]. (A similar mechanism is 
operative in the 5-dimensional composite Higgs models [7], where the role of the top 
partner is played by the Kaluza-Klein excitations of the top.) In this paper, we will 
focus on a minimal model that incorporates the top Yukawa via collective symmetry 
breaking and explicitly realizes the structure of Eq. (1.1). We will present direct and 
indirect bounds on the model and discuss their implications for naturalness in light of 
the 125 GeV Higgs. We will also consider predictions for the deviations of the Higgs 
and top properties from the SM. 

Our model is basically identical to the top sector of the Littlest Higgs [5], and 
we will make use of many results derived in the context of that model. The original 
Littlest Higgs is severely constrained by precision electroweak data [8]. The constraints 
come almost entirely from the extra gauge bosons of the model, whose masses are 
required to be above 2-3 TeV. In itself, this is not a problem for naturalness. However, 
the structure of the Littlest Higgs imposes a tight relation between the gauge boson 
and top partner masses, so that multi-TeV top partners are required, which in turn 
implies strong fine-tuning. This problem can be avoided by modifying the model, 
by introducing an additional symmetry (T-parity) to forbid tree-level corrections to 
precision electroweak observables [9], by decoupling the top and gauge boson partner 
mass scales [10], or simply by slightly lowering the cutoff and getting rid of the extra 
gauge bosons altogether [11]. Thus, while the structure of the top sector is robust 
- it is in effect fixed by the naturalness requirement - the gauge and scalar sectors 
appear quite mo del- dependent, both in their structure and in the associated mass 
scale. Motivated by these considerations, we consider the top sector in isolation, and 
identify the predictions that are in a sense unavoidable once the cancellation mechanism 
in Fig. 1 is postulated. This approach is similar to the bottom-up attitude to SUSY 
phenomenology advocated in Refs. [2]. 

The rest of the paper is organized as follows. The minimal model for the fermionic 
top partner is presented in Section 2. Section 3 discusses naturalness of electroweak 
symmetry breaking in this model, assuming a 125 GeV Higgs boson. Section 4 summa- 
rizes existing experimental constraints on the model, divided in three groups: precision 
electroweak, flavor constraints, and direct searches at the LHC. Sections 5 and 6 dis- 
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cuss the expected deviations of the Higgs and top properties, respectively, from the SM 
predictions. We summarize our findings and conclude in Section 7. A number of useful 
formulas are collected in the Appendix. 



2 Minimal Model for Fermionic Top Partner 

We begin with a non-linear sigma model describing spontaneous SU (3) — >■ SU (2) global 
symmetry breaking by a fundamental vev. The sigma field is 



V = «p ^ , (2.1) 




where t a are the broken generators (a = 1 ... 5), n a are the corresponding Goldstone 
bosons, and / is the symmetry breaking scale (we assume f & 1 TeV). We identify 
the £77(2) doublet of Goldstone bosons with the SM Higgs doublet H, and ignore the 
remaining one which plays no role in our analysis: 

" A =U o)' < 2 - 2 ' 

To generate a top Yukawa coupling without introducing one-loop quadratic divergences, 
we introduce an SU(3) triplet of left-handed Weyl fermions, \l = (& 2 Q, U)\, and two 
£77(3) singlet right-handed Weyl fermions, ur and Ur. Here Ql = These 
fields are coupled via [5, 12] 

£ = -\ iUr V^ X l - MfU R U L + h.c. (2.3) 

Expanding the sigma field up to terms of order 1/f 2 gives 

C = -f (Xiu R + \ 2 U B ) ] U l - X iUr HQ l + j-(rfH)u R U L + h.c. + . . . (2.4) 
where H = {ia^Hy. The fermion mass eigenstates are 

m TT rr >>lUR + \ 2 U R 

J-l — Ul, J-r 



2 



, X2UR - XiUr . , 

tL = UL ' tR= ^x{ + x> ' (2 ' 5) 



where we neglected the Higgs vev v, assuming v <g /. (The mixing angles and masses 
with full v dependence are given in Appendix A.) We identify t = (t L , t fl ) with the SM 
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top quark, b L with the SM left-handed bottom, and T = (T L , Tr) with the top partner, 
whose mass is 



m T = ^\\ + \lf. (2.6) 
The interaction terms in the mass eigenbasis become 



Ant = -XA&Ql ~ \tT ] b HQ l + ^-(H^H)T R T L + ^(H^H)t R T L + h.c. + 



(2.7) 



where we defined 



V + A 2 v Ai + A 2 



The first term is simply the SM top Yukawa; the next two terms reproduce Eq. (1.1), 
ensuring the cancellation of the one-loop quadratic divergence (note that A^ = A^ + A^); 
while the last term does not contribute to the Higgs mass renormalization at one loop, 
and thus does not spoil the cancellation. The cancellation is also easy to understand in 
terms of symmetries of the model: the first term in (2.3) preserves the global SU(3), 
so that in the limit X 2 — > the Higgs is an exact Goldstone boson and is therefore 
massless. On the other hand, the second term in (2.3) breaks the SU(3) explicitly, 
but it does not involve the Higgs at all, and so cannot generate the Higgs mass on its 
own. Thus, both couplings need to enter any diagram contributing to the Higgs mass 
renormalization, and at the one-loop level the diagrams involving both A's are at most 
logarithmically divergent. 

The Higgs can be given its usual SM gauge couplings by weakly gauging the SU (2) x 
U(l) subgroup of the SU(3). As explained in the Introduction, we do not consider 
extended gauge sectors here: the gauge structure of our model is the same as SM. 
The new top-sector fields Ul, Ur have the same gauge quantum numbers as the SM 
right-handed top, (3, 1)4/3. 

Non-linear sigma model interactions become strongly coupled at a scale A as 47r/, 
where another layer of new physics must occur. The effects of that physics on weak- 
scale observables can be parametrized by adding operators of mass dimension > 4, 
suppressed by appropriate powers of A, to the lagrangian. The leading (dimension-6) 
operators are 

£uv = j^jy 2 {V^D.V) 2 + ^W^(V^V) , (2.9) 

where is the covariant derivative including the SU(2) x U(l) gauge fields; W and B 
are the SU(2) and U(l) field strength tensors, respectively; c\ and c 2 are dimensionless 
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coefficients, which are unknown but expected to be of order 1; and 

Q a = ( ° a ) • (2.10) 

The two operators in Eq. (2.9) contribute to the T and S parameters, respectively, in 
precision electroweak fits (see Sec. 4.1). We do not include operators involving the top 
quark, since they are not strongly constrained at present. 



3 Higgs Mass and Naturalness 

An appealing feature of the class of models we're dealing with is a simple, rather 
predictive description of the electroweak symmetry breaking (EWSB). At tree level, 
the Higgs is a Goldstone boson and the Higgs mass parameter /i 2 = 0. At one loop, the 
leading (log-divergent) contribution to the Higgs mass parameter from the diagrams in 
Fig. 1 is given by 

\ 2 m 2 A 2 

6^ = -3^f\o g \. (3.1) 

Naive dimensional analysis (NDA) suggests that this is the dominant contribution to the 
Higgs mass: two-loop quadratically divergent contributions are suppressed by a power 
of log(47r) 2 w 5, while gauge boson loops (assuming that their quadratic divergences 
are canceled at a scale close to 1 TeV) are down by (g/X) 2 . Note that Eq. (3.1) 
automatically has the right (negative) sign to trigger EWSB. 

If the LHC hint is correct and there is indeed a 125 GeV Higgs boson, then /z 2 can 
be treated as known, since = y/2\fj,\. In our model, this essentially fixes the top 
partner mass, up to logarithmic dependence on A. For definiteness, we take A = 4nf; 
to leading order in v/ f, 

2tt , , 

A « — m T sin2a, (3.2) 
At 

where a is the mixing angle in the right-handed top sector (at leading order at v/f, 
tana = Ai/A 2 ). For example, for a = 7r/4, we obtain 

m T « 236 GeV. (3.3) 

Unfortunately, the top partner at this mass is excluded by precision electroweak con- 
straints, see Section 4. The mild a dependence does not change this conclusion if a is 
varied within a reasonable range. 

The only way to raise and salvage the model is to assume that the gauge-loop 
and/or two- loop contributions to fi 2 are enhanced, and partially cancel the top- loop 
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Figure 2. Fine-tuning as a function of the top partner mass (in GeV). We fixed m;, = 125 
GeV. The band corresponds to varying the mixing angle a between 0.2 and 1.1. 



contribution. 1 This requires a certain degree of fine-tuning; we quantify it by defining 

A = M (3.4) 

^obs 

where /i Q b s = rn h /\/2 w 88 GeV. Required fine-tuning as a function of the top partner 
mass is shown in Fig. 2, where the band corresponds to varying the mixing angle 
a between 0.2 and 1.1, corresponding roughly to the range where both Ai and A2 
are perturbative. This plot should be kept in mind as we discuss the experimental 
constraints on the model below. 

A Higgs quartic coupling, A ~ 0.12, is required to accommodate the Higgs vev v = 
246 GeV along with a 125 GeV mass. In our model, there is no tree-level quartic, but 
at one loop the quartic is generated by quadratically divergent terms in the Coleman- 
Weinberg potential [5, 14]. In our minimal model, the quartic generated by the top- 
sector is in fact only logarithmically sensitive to the cutoff, and is thus expected to 
be small. However, the contributions to global symmetry breaking due to gauging the 
SM SU{2)l x U{1)y do generate quadratically divergent contributions to the quartic. 
These diagrams are dominated by physics at the scale A, and hence cannot be computed 



1 In fact, Ref. [13] argued that the two-loop contribution in the Littlest Higgs is enhanced compared 
to the NDA estimate, and estimated that it is of the same order as the logarithmically divergent one- 
loop contribution. Since the two-loop contribution is UV-dominated, its magnitude (and sign) cannot 
be determined without specifying a UV completion and performing a calculation in a UV-complete 
model. 
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without specifying a UV completion, but NDA estimates show that an 0(g 2 ) quartic 
can be generated without tuning. 



4 Experimental Constraints 

The model in Eq. (2.3) has three parameters: the symmetry breaking scale / and 
two dimensionless couplings Ai,2- One combination of the couplings has to be fixed 
to reproduce the known top Yukawa, leaving two independent parameters. In our 
discussion of experimental constraints, we will use the top partner mass vht and the 
rotation angle a between the gauge and mass eigenstates in the right-handed fermion 
sector. That is, a is defined by 



tji = cos a Ur — sin a ur , T R = sin a% + cos aU^. 



(4.1) 



The relation between (%, a) and the Lagrangian parameters, to leading order inv/f, 
is given in Eqs. (2.5), (2.6). In the analysis below, we will use generalizations of these 
formulas to all orders in v/f, see Appendix A. It is also worth noting that at order 
v/f, mixing between the left-handed fermion fields ul and Ul is induced; the mixing 
angle (5 is approximately given by 



sin /3 « tan a 



rriT 



(4.2) 



Again, we will use the exact expression for this mixing angle, given in Appendix A. 
This mixing induces the off-diagonal vector boson couplings to fermions, ZtT and 
WbT, which play a crucial role in the phenomenology of the model. Both couplings 
are proportional to sin f3. 

4.1 Precision Electroweak Constraints 

The top partner T does not induce tree-level contributions to precision electroweak ob- 
servables. At one-loop, oblique corrections to the electroweak gauge boson propagators 
induced by diagrams involving the T are given by [15] 



St 
T t 
U t 



2tt 
3 



3^ 



rn 



16n sljClj rn 



-.2 r 



2tt 



logo;* + cl 



[l + x t 



r s 2 

JL 



-l-ci- 



1 - x t 



+ 



2clxf (3 - x t )\ogx t 



:i-xt) 



8c 2 1 
~3 



log x t 



Sp logxt + c 



2 (1 + X t ) 

Pf l-x t ) 



2cpXt(3 - x t )logx t 
(l-x t f 



H 1 



(4.3) 
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where x t = m 2 /m^, and s w is the sine of the Weinberg angle. In addition, there is a 
contribution due to the shift of the Higgs couplings to the electroweak gauge bosons 
from their SM values [16]: 



s h = 


1 


log 


m h 




A 


T h = 


3 m 2 w 


lo£ 


m h 




' A 



(4.4) 

where A is the scale where the Higgs loops are cut off. We will assume A = Airf. 
Furthermore, the operators induced by the new physics at scale A, given in Eq. (2.9), 
contribute [17] 



S 



uv 



^c s m 2 w 

2 f2 



ng 2 f 



ctm 2 w 



Tvv = - v ^j- 2 . (4.5) 

The only important non-flavor-universal correction is the top-partner loop contribution 
to the ZbJ)L vertex. To leading order in the limit rrix ^> m t ^> mw, this is given by [15] 

. h z g a mf A?, mt nS 
c w 8nsi m z w m z T Y 2 mi 

The correction to the Zh^bn vertex is negligible since it is not enhanced by the top 
Yukawa coupling. 

The results of a fit to the precision electroweak observables [18] are shown in 
Fig. 3, where we also included contours of constant fine-tuning computed according 
to Eq. (3.4). We conclude that: 

• The lower bound on the top partner mass from precision electroweak observables 
is approximately 500 GeV; 

• The corresponding minimum level of fine-tuning on the Higgs mass is about 20%. 
This is significantly better than in the MSSM with a 125 GeV Higgs, and com- 
parable to the NMSSM with large A [3]; 

• These conclusions do not depend strongly on the operators induced by the UV 
completion of the model, as long as the size of these operators is roughly consistent 
with naive dimensional analysis. 
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Figure 3. Precision electroweak constraints on the minimal fermionic top partner model, in 
the (rriT,Oi) plane. The three panels display the variation of the bounds depending on the 
coefficients of the UV operators: (a) c s = a = 0; (b) c s = ct = 1; (c) c s = +1, ct = —1. Thick 
black lines represent constant fine-tuning contours: from left to right, 20%, 5%, and 2% fine 
tuning. 



4.2 Flavor Constraints 

By selecting the top quark to be the only one with a partner, and introducing mixing 
between the SM top and its partner, our model explicitly breaks the approximate flavor 
symmetry of the SM, leading to potential constraints from flavor-changing processes. 
We investigate these constraints in this section. 

Including the mixing between the three SM generations, the mass terms of the 
up- type quarks in the gauge basis form a 4 x 4 matrix M^, while the down-type mass 
terms are described by a 3 x 3 matrix M% . (Here and below, capital indices run from 
1 to 4, and the lower case indices from 1 to 3.) Diagonalizing these matrices requires 

M u -> L U M U R\ , 

Md — >■ L d M d R\ , (4.7) 

where L and R matrices rotate the left-handed and right-handed quark fields, respec- 
tively. The charged-current interactions in the gauge basis have the form 

C c .c.=gW-J + »+ c.c, 

J + » = -^UinPsYiiD^ , (4.8) 
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where 



In the mass basis, the charged current becomes 

# = ^M4)/(p 3 )$ML«* , (4-io) 

so that the generalization of the CKM matrix in our model is 

(Vckm)] = (Li)f(P 3 ) k K (L d ){. (4.11) 

The elements of this matrix should in principle be determined by a fit to data. We 
will not attempt such a fit here. Since the SM CKM matrix provides an excellent 
description of flavor-changing processes for the first two generations and the b quark, 
we assume the following structure: 



V, 



CKM 



( V u d V us V u b \ 
V c d V cs V c b 
cpV td cpVu cpVtb 
\ -SfjVu -spV ts -s p V tb J 



(4.12) 



where are SM CKM elements. With this assumption, all flavor- violating new physics 
effects in K and B systems appear at loop-level only. 

Unlike the SM, rotations (4.7) induce tree-level flavor- changing neutral currents 
(FCNC) in the left-handed sector [19], since the weak-singlet Ul mixes with the SM 
up-type quarks. The Z boson couples to the current 

^ NC = (UlYY (T 3 - slQ) 3 ; (U l )j, (4.13) 

where 



Rotation to the mass basis yields flavor-changing couplings, proportional to 

V FCNC = Ll(T 3 -s 2 w Q)L u . (4.15) 

These can generate tree-level contributions to rare D meson decays and anomalous 
D — D mixing, and flavor-changing top decays. Such contributions are however 
completely absent if 

12x2 

L„ = | c p -s p | , ( LU>) 

S/3 Cp 
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since the only flavor- violating Z coupling in this case is ZtT. Eq. (4.11) then requires 
Ld = Vckm- We will assume this texture in our analysis. Note, however, that due to 
large theoretical uncertainties associated with the D system and the highly suppressed 
rates for anomalous top decays, significant deviations from this texture can still be 
consistent with experimental constraints [19]. 

At the one-loop level, our model predicts new contributions to AF = 2 and AF = 1 
processes in B and K systems. Let us first consider AF = 2. The effective Hamiltonian 
that governs the B° s — system is 

H Bs = y^ M w Yl h\jF{x I) x,j-M w )x(bs) (v ^ A) {sb) {v „ A) , (4.17) 

I,J=u,c,t,T 

where we defined A/ = Vf b Vi s and xi = Mj/Myy- The F functions are given in 
the Appendix B. Hamiltonians for the K° and B^ systems are obtained by simple 
substitutions, b — > d and s — > d, respectively. At leading order in v/f expansion, 
our results agree with Refs. [20-22]; however, our expressions are exact in v/f. To a 
good approximation, the size of the new physics effects in AF = 2 observables can be 
estimated as the fractional deviation of the Wilson coefficient in Eq. (4.17) from its SM 
value. (This estimate does not take into account some effects, such as the running of 
the Wilson coefficient between the scales mx and m t , which are however expected to 
be small.) We find that the maximum deviations on the parameter space of our model 
are: 0.5% for Am^; about 20% for ex] and about 35% for \Am(Bd)\ and \Am(B s )\. 
Such deviations are currently easily allowed by data: see, for example, [23]. 

We next consider the two most constrained AF = 1 decays, b — » 37 and B s — > 
H + H~ . The b — > 37 amplitude is proportional, in the leading- log approximation, to the 
coefficient C 7 of the operator P 7 = ^2m h {sLa^ v bR)F^ v1 evaluated at the scale m b . The 
top-quark contribution to this coefficient is given by 

Xt = -\Al(x t ) V ^ - \Fl{x t ) (V 4/23 - , (4-18) 

where the functions A^x) and F^(x) can be found in Appendix B, and r] = a s {m t ) / a s {rrib) 
The only effect of the top partner is to replace 

Al{x t ) -> cIAKx^ + s^AKxt), 

F$(x t ) -> c 2 ^(x t ) + *j,F*(xr), (4.19) 

in these expressions. (The first term in the v/f expansion of these formulas agrees 
with Refs. [22, 24]; however, our formulas are exact in v/f.) The resulting deviations 
of the b — > 57 branching ratio from the SM are shown in the left panel of Fig. 4. 
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Figure 4. Fractional deviations of the B — > X s ^f (left panel) and B s — > /x + pT (right panel) 
branching ratios from the SM predictions. The thick line on the right panel corresponds to 
the LHCb upper bound on Br(S<j — > fj, + /j,~); all points above the line are now ruled out. 



In the region of interest, the deviations are at most about 5%. Given that both the 
experimental measurement [18] and the NNLO SM theoretical prediction [25] have 
uncertainties between 5 and 10%, such deviations cannot be currently ruled out. The 
right panel of the figure shows the deviation of the B s — > fi + fi~ branching ratio from 
the SM prediction, evaluated using the formulas given in Ref. [26]. We also indicate 
the region ruled out by the recent LHCb bound [27], Bi(B s — > < 4.5 x 10 -9 at 

95% c.l., which is only a factor of 1.5 above the SM prediction. This is the strongest 
current bound on the top partner from flavor physics, even though it is still weaker 
than precision electroweak constraints. Notice that the results of Ref. [26] are valid to 
leading order in the v/f expansion. Given the potential importance of this bound, a 
more precise calculation is desirable. 

4.3 Direct Searches at the LHC 

The two production mechanisms for the top partner are strong pair-production, qq/gg — > 
TT, and electroweak single production, bq — > Tq' or qq' — > Tb. The production cross 
sections at the 7 TeV LHC are shown in Fig. 5. For pair-production the cross-section 
is calculated at NNLO in QCD using Hathor vl .2 [28], with renormalization and fac- 
torization scales set to the top partner mass. For single-production the cross-section is 
calculated at LO using MadGraph5 vl.3.32 [29]. While an NLO calculation of single- 
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Figure 5. Production cross section of top partners at the 7 TeV LHC. For pair-production 
(red), we use Hathor vl . 2 [28] to calculate at NNLO in QCD. For single production (blue), we 
use MadGraph5 vl.3.32 [29] to calculate at LO. The single production cross-section depends 
on a, with the band indicating the cross-sections for 0.2 < a < 1.1. 



top production cross section is available [30], its use is not justified in out study due to 
large model uncertainty in the bTW coupling. In this case, we use the MadGraph default 
setting for renormalization and factorization scale, variable event-by-event. (For both 
pair- and single-production, varying renormalization and factorization scales within a 
factor of 2 leads to at most a few % variations in the cross sections.) At the 7 TeV LHC, 
due to the relatively small phase space for producing heavy particles, single production 
overcomes its electroweak suppression and can be comparable to pair production. 

Decay channels of the top partner include th, tZ and bW [12, 31]. In the limit 
/ S> v, the branching ratios are 25%, 25%, and 50%, respectively, as can be easily seen 
from the Goldstone boson equivalence theorem. An explicit calculation of the partial 
widths yields [12]: 



m) = Trf- / (**, x h ) [(1 + x\ - x\) ((C' th ) 2 + (C£J 2 ) + ^ th C* h x t ] , (4.20) 
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Final state Raw Br^A, (%) Rescaled Br^A (%) BrjBr^ejAj (%) 



bWbW 
bWtZ 
bWtH 

tztz 

tHtH 



0.36 
0.034 
0.022 
0.0015 
9.9 x 10~ 4 



0.29 
0.027 
0.018 
0.0012 

7.9 x 10~ 4 



0.12 
0.0046 
0.0011 

8.5 x 10~ 5 
7.5 x lO" 6 



Table 1. Estimated raw Br^ejAj, rescaled Br^ejAj and Br^Brjej^j for the various decays of 
pair produced TT. All values assume = 400 GeV and a = 7r/4. See text for the definition 
of raw and rescaled efficiencies. Bere Br^ denotes the dileptonic branching fraction for WW, 
which is common to all decay modes. 

where X{ = rni/rriT, the kinematic functions are defined as 



and the constants appearing in the tTh vertex are given in Appendix A. 

There exist several searches for vector-like top partners at CMS and ATLAS [32- 
37]. These searches focus on pair production and on one particular decay mode of the 
top partner, either bW or tZ, and assume 100% branching fraction to that mode. In 
our model, the signal is generally a mixture of pair and single production, and multiple 
decay channels are possible. As a result, the bounds on the top partner masses obtained 
by CMS and ATLAS are not directly applicable, but it is possible to "recast" the 
published analyses to estimate the bounds in our model. 2 Below we present such an 
estimate, based on the CMS search in the bbWW final state with 5.0 fb _1 integrated 
luminosity [32]. In the interesting parameter space of our model, the dominant decay 
mode for the T is bW, making bbWW searches most sensitive. Furthermore, the CMS 
analysis places the strongest bounds as it is updated to use the full 2011 dataset. 

The number of signal events expected in a given analysis can be written as 



where <7j is the cross-section for each production channel, C is the integrated luminosity 
used in the search, Br^ is the branching fraction for an event produced via channel i 
to result in the final state j after the decay of all unstable particles, is the efficiency 

2 See [38] for a recent theoretical analysis that recasts experimental searches in terms of limits on 
top-partners with general values of the branching fractions to bW, tZ, and th. 
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for detecting the final state j in a given analysis, and A^ is the acceptance for the final 
state j. (Note that e and A depend on the production channel, since final-state par- 
ticles have different kinematic distributions depending on the production mechanism.) 
The efficiency and acceptance for the particular production and decay mode assumed 
in the CMS analysis (TT, with T — > bW) can be found in [32]. We estimated the eA 
of all other relevant final states by modeling the acceptance and selection cuts of [32] 
on a sample of Monte Carlo (MC)-generated TT events. (Since the CMS analysis re- 
quired two isolated leptons, and vetoed events with dilepton invariant mass close to 
the Z boson, the efficiencies for events with a single T to pass the cuts are extremely 
small, and we did not include the single production channel in our analysis.) For this 
estimate, we generated parton- level events using MadGraph 5 vl.3.32 [29], showered 
and hadronized them using Pythia 6.426 [39], and applied simplified detector sim- 
ulation using PGS 4 . [40] . Unfortunately, PGS 4 . significantly underestimates the 
efficiency of b-tagging, compared to the TCHEM algorithm used by the CMS in this 
analysis. 3 To address this issue, we ignored the b-tag information provided by PGS, 
and instead applied pr-dependent TCHEM efficiencies [41] to the b-jets in our sample. 
This procedure yields the "raw" eA values for all possible final states, as a function of 
rriT and a. For example, values of eA for inx = 400 GeV, a = tt/A are listed in the 
first column of Table 1. The MC simulation and analysis procedure was validated on a 
sample of events with the final state considered by CMS, WbWb with 2 leptonic Ws. 
We found that the eA determined using our simulation for a 400 GeV top partner is 
0.36%, compared with the value of 0.29% quoted in the CMS analysis. Given the crude 
nature of our MC simulations, this level of agreement is very reasonable. Even so, in 
deriving the bounds, we rescale the raw MC estimates of BieA by a correction factor 
of 0.29/0.36; in other words, we use the eA quoted by CMS for the WbWb channel, 
and use the MC to estimate the relative eA of other channels with respect to WbWb. 
The resulting estimates are collected in Table 1. It is clear that the rates of events 
from final states other than WbWb that pass the analysis cuts are quite small. While 
our estimates of those rates suffer from significant systematic uncertainties due to the 
crude detector simulation used, it is reassuring that even if the rates were inflated by 
a factor of two they would remain subdominant. Thus, our bounds on the top partner 
mass are robust. 

The estimated 95% C.L. exclusion region as a function of my and a is presented in 
Fig. 6. The typical bound on the top partner mass is about 450 GeV, with somewhat 
weaker bounds for small a. 4 It is clear that direct collider searches are just beginning 

3 This can be easily seen by comparing PGS and TCHEM efficiencies on an SM bb sample. The 
peak efficiencies are 0.4 for PGS and 0.7 for TCHEM. 

4 A reanalysis of the published LHC searches in the context of the "Bestest" Little Higgs model 
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is 




200 300 400 500 600 

m T 

Figure 6. The estimated 95% excluded region from the CMS T — > Wb search [32] in terms 
of Mt and a is shown as the red contour. The green is the allowed region. The thick black 
line corresponds to 20% fine-tuning. The CMS analysis does not quote efficiencies below top 
partner masses of 350 GeV, so no bound is shown (grey region). This low-mx region is in 
any case already ruled out by precision electroweak fits. 

to probe the region that is not already ruled out by precision electroweak constraints. 
Note that the least fine-tuned parameter space regions will be probed by direct LHC 
searches in 2012. 

ATLAS has also published a search for a singly-produced vector-like heavy quark, 
in particular using the channel pp — > Qq — > Wqq' [43]. The lower bound on the Q mass 
of order 900 GeV was reported. This search potentially has sensitivity to our model, 
since single top partner production contributes to this final state. However, the Q in the 
ATLAS analysis was assumed to have direct coupling to first-generation quarks, while 
in our model T only couples to the third generation, resulting in much lower production 
cross sections. As a result, this analysis does not yet put interesting bounds on the 
top partner mass. For example, ATLAS sets a bound on a{pp — > Qq)xBi(Q — > Wq) 
of about 2 pb for a 500 GeV Q mass; in our model, a single top-partner production 
cross section of this size only occurs at the upper edge of the band shown in Fig. 5. 
This requires values of a ~ 1, which are already ruled out by precision electroweak 

appeared recently in Ref. [42], where similar bounds on the top partner mass were found. 
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constraints. Note, however, that the published ATLAS analysis only uses 1.04 fb _1 of 
data; updated versions of this analysis will have sensitivity to interesting regions of the 
top partner parameter space. 

5 Higgs Properties 

If the LHC evidence for the Higgs boson at 125 GeV is correct, detailed measurements 
of the Higgs production cross section and branching ratios should be possible within the 
next few years. In our model, these properties deviate from the SM predictions. There 
are two important effects. First, the hWW and hZZ couplings are shifted, 5 leading to 
deviations in the WW and ZZ branching fractions and, via the WMoop contribution, 
in the Br(/j — > 77). Second, loops of top partners produce corrections to the hgg and 
fo<yry vertices, leading to deviations in the expected production cross section and, again, 
Br(h -)■ 77). 

The production rate of h via gluon fusion is proportional to T(h — > gg). Assuming 
that gluon fusion is the dominant Higgs production mechanism, the rates a(pp — > h — > 
VV) in our model, normalized to their SM values, are 

T(h -> gg)Br(h -» VV) 

V T SM (h -> gg)Br SM (h -+VV)' 1 ' ] 

where V = 7, Z, W. The total Higgs decay rate at = 125 GeV is dominated by 
the bb mode. The bottom Yukawa coupling can be incorporated in our model as an 
explicit breaking of the global symmetry; this would not spoil naturalness due to the 
small numerical value of i/b- At leading order in v/f, this results in the hbb coupling 
identical to the SM value. There may be corrections at higher orders in v/f; however, 
their form is not fixed by the symmetry, and is model-dependent. If they are ignored, 
we simply get 

T(h -» gg)T(h -> VV) 

V T SM (h^gg)T SM (h^VV)- { ' ] 

Note that the dropped terms in the hbb vertex are potentially of the same order as the 
corrections to the hgg and hVV couplings, so these predictions have an inherent 0(1) 
ambiguity. Still, we compute them as an indication of the likely size of the effect. We 
should also note that our predictions for ratios of rates, such as for example R^/Rw, 
are free of this ambiguity. 



5 These shifts are due simply to the composite nature of the Higgs, and not to the presence of the 
top partners. They can be described within the framework developed in Refs. [44]. 
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The ratios of the h — > WW and h — > ZZ decay rates to the SM predictions are 
given by 

r(h -> w + w~) r(h -> zz) 



r SM (/i -> w+w- 

The /i — )• 77 decay rate is [45, 46] 



r SM (fr -> zz) 



2m|, 
9 2 P 



(5.3) 



r(/i 77) 



mi 



1024tt 3 



fiW + ^Q-A^Fv^r,) + 3Q^U 



top 



(5.4) 



where Tj = 4m 2 /m\; the sum runs over all SM fermions except the top; Qi is the electric 
charge of the i-th fermion and N c ^ its color multiplicity (3 for quarks, 1 for leptons). 
The top contribution to the decay amplitude in our model is given by 



A 



V2m 



w 



top 



gm t 



CtthFl/2{Tt) H -CTThFl/2{T T ) ) 

mT J 



where the constants C are given in Appendix A; while in the SM, .Atop 
we used the standard notation for the loop functions, 

= 2 + 3x + 3x(2 - x)f{x) , 
F 1/2 (x) = -2x(l + (l-x)f(x)) , 

'1 



(5.5) 

Fi/ 2 (r t ). Here 



/(*) 



sin 



if x > 1, 




Z7T 



if x < 1. 



(5.6) 



The h gg decay rate is given by [45] 



m 3 h 

5127T 3 myy 



top 



(5.7) 



The predicted rates Rw = Rz and i? 7 are shown in Fig. 7. Comparing with the 
precision electroweak constraints, we conclude that large suppression of the rates in 
both WW/ZZ and 77 channels is possible: the rates can be as low as 30% of the SM 
prediction. Deviations are the strongest for the least fine-tuned regions of parameter 
space: for example, if we demand EWSB fine-tuning of 5% or better, the minimal 
possible deviation in Rw and R 1 is 20%. As noted above, these predictions should be 
taken with a grain of salt, since they can be modified by the model-dependent 0(v/f) 
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terms in the hbb coupling. Still, it is interesting that large, potentially observable 
deviations from the SM may occur throughout the natural parameter space. 

As remarked above, the ratio R^/R\y provides a robust test of the structure since 
it's insensitive to the model-dependent embedding of the bottom Yukawa. Unfortu- 
nately, throughout the parameter space of our model, the deviations of this ratio from 
the SM prediction are well below 1%, too small to be observed. The reason is that to a 
very good approximation, the fractional deviations of the hWW coupling and the top 
loop contributions to hgg and h'yy are the same. 

6 Top Properties 

At order v/f, the lighter top eigenstate, which we identified with the SM top, actually 
contains an admixture of the S77(2)-singlet left-handed field Ul- As a result, the chiral 
structure of the top couplings to the Z deviates from the SM predictions at this order. 
To quantify this effect, in Fig. 8 we plot the ratio of the vector and axial components 
of the tiZ coupling expected in our model, normalized to their SM values. Deviations 
of order 10% or more in g^, and up to 30% in gy, are possible in regions consistent 
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with precision electroweak constraints. It was estimated that the 14 TeV LHC with 
3000 fb _1 integrated luminosity would be able to probe g& at the 5-10% level and gy 
at the 15-30% level [47]. A proposed 500 GeV linear electron-positron collider would 
reach precision of 2% on g^ and 5% on g v [48, 49]. Though we would expect the top 
partner to be discovered in direct searches before these measurement become possible, 
they would still be of great interest to confirm the structure of the model. 

7 Conclusions 

Naturalness, together with evidence that electroweak-symmetry breaking sector re- 
mains weakly coupled up to scales well above 1 TeV, implies that a light Higgs must 
be accompanied by new particles that cancel the quadratically divergent Higgs mass 
contribution from the SM top loop. In SUSY, these particles are scalar top (stop) 
quarks, and LHC phenomenology of stops has been a subject of much work recently. 
In this paper, we studied an alternative which has not received as much attention so 
far: naturalness restoration by spin-1/2 top partners. We focused on a minimal model 
where this mechanism is realized, which is essentially the top sector of the Littlest Higgs 
model. We explored current experimental constraints on this model from all relevant 
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sources: precision electroweak fits, flavor physics, and direct LHC searches. We found 
that the current bound on the top partner mass is about 500 GeV, and is dominated by 
precision electroweak data, although direct searches are rapidly entering the hitherto 
allowed mass range. Given these bounds, accommodating a 125 GeV Higgs boson in 
this model requires only a modest level of fine-tuning, of order 20%. Thus, we conclude 
that natural EWSB is possible in theories with sub-TeV-scale spin-1/2 top partners. 

In the near future, direct searches for the top partners at the LHC will continue, 
gaining more sensitivity as more data is collected. The decay channels of the top 
partner include bW, tZ and th, all of which have order-one branching ratios; this 
situation is not special to our model but is in fact quite generic. Also, while existing 
searches focus on pair-production of the top partners, in our model single production 
dominates in parts of the parameter space. To maximize sensitivity to top partners, 
experiments should extend the menu of searches to encompass all available production 
and decay modes. Another interesting handle not used in the top partner searches so 
far is jet substructure: the top partner decay products, such as t, Z, W and h, are 
typically relativistic in the lab frame in the relevant mass range, so that their hadronic 
decays can be identified as jets with unusual substructure. Recent phenomenological 
studies [50, 51] show interesting potential of such searches. 

As a complementary handle, measurements of the Higgs and top properties at the 
LHC may be sensitive to deviations from the SM predicted by our model. While these 
predictions are quite model-dependent, our study indicates that large deviations in 
h WW/ZZ and h — > 77 rates are possible. 
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A Masses, Mixing Angles and Couplings of the Top and Its 
Partner 



Ignoring the Goldstone fields that are eaten by the SM gauge bosons after EWSB, the 
sigma field V has the form 

V = e W (™I)( ° f ) , (A.l) 

where 

' 1 
, 1 

and 

1 v + h 



(A.2) 



(A.3) 



V2 f 

Here v is the Higgs vev, and h is the physical Higgs boson. The exponent can be easily 
expanded using the fact that I 2 = 1: 

exp(zaf) = cos a + i sin a / . (A. 4) 

The kinetic term of the sigma model has the form 

A™ = (D^V) ] (D^V) , (A.5) 

where is the covariant derivative. This term contains a canonically normalized 
kinetic term for the Higgs, as well as masses for the SM gauge bosons; in particular, 

m w = \d 2 f 2 sin2 s > ( A - 6 ) 

where we defined a = v / (a/2/). The measured value of m w can be used to compute v 
from this formula; in the limit / — > oo, v tends to its SM value, 246 GeV. 
Using Eq. (A. 4), the top mass terms take the form 

Anass = (4 U ] R )M ( ^ ) + h.C. , (A.7) 



L 



where 



M = f( AlSin " Al ^° S " ). (A.8) 
A2 



Diagonalizing Ai^M, we find the masses of the top quark t and its partner T: 



, (Af + \l)P I I 4AfA|sin 2 a 
m UT = - n I l=b (Aj + A|)2 1 • (A.9) 
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The rotation between gauge eigenstates (u, U) and mass eigenstates (t, T) is given by 



ti = cos (3 ul — sin/3 Ul, T l = sin (3 ul + cos (3 Ul ', 

in = cos auR — sin aUn, Tr — sin ait,R + cos a Ur , (A. 10) 



and the mixing angles are 



1 i 2AiA 2 cos a 
a — - tan 



2 A 2 — A 2 

1 i A? sin 2d 
— tan 1 

2 A 2 , + A 2 cos 2a' 



P= otan- 1 -^-^— — . (A.11) 



Mass and mixing angle formulas quoted in the main text are obtained by expanding in 
v/f and keeping the leading order terms only. 

It is also useful to invert these formulas and express the Lagrangian parameters 
(Ai, A2, /) in terms of physical parameters (m t , mr, a): 

( \/2mw \ l/ 2 , • 2 \i/2 / . 2 2 \V2 

j = I I (cos a + x t sin a J (sin a + x t cos a) 



x] /2 



Ai 



/ gmt A 1 

\\/2m w J (cos 2 a + x i sin 2 a) 1 



/2 



^=(-^"1 1 T75, (A.12) 

\y/2m w J (sin 2 a + x t cos 2 a) 7 

where = m 2 /m^. For example, together with the second line of Eq. (A. 11), this 
expressions give the angle (3 in terms of the physical parameters, which was used in the 
calculation of precision electroweak parameters in Sec. 4.1: 

1/2 

sin/3 = -^^=. (A. 13) 

V cot a + x t 

The couplings of the top and its partner to electroweak gauge bosons are given by 

£ g = — t— b\o* (cos/3t L + sin/3T L ) W~ + c.c. 
V2s w 

+ ((/ tt£ t\oH L + a ttR + 9ttl Tla"T L + T R a^T R ^j 

+ (? tTL ^T L ^+ c.c, (A. 14) 
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where 



e / cos 2 /3 2s 2 \ 2es 



9ttL — — n 7T ; 9ttR 



i w 



c s \ 2 1 / ' y V ' 

e / sin 2 (3 2s 2 \ 2es^ 



9ttl — I — ~ 77- I ; 9ttr „ 

^tl = ^— • (A. 15) 

Their couplings to the Higgs boson are 

£ yu k = ~ (Cttht ] L t R + CtthTiTr + C^ th t ] R T L + C^TjjtL^j -j= + h.c, (A. 16) 
where 



CWi = Aj cos a cos(a — CVt/i = — Ai sin a sin(a — /?); 
Crt/i = — Ai cos a sin(a — (5)\ C^ th = Ai sina cos(a — /3). (A. 17) 

Finally, the Higgs boson coupling to the electroweak gauge bosons are given by 

C h w = 2 cos a (m 2 w W + »W- + \m 2 z Z»Z^\ ^- . (A. 18) 

\ z J Zmw 

These couplings are suppressed compared to the SM values by a common factor, 



(A.19) 

B Loop Functions Appearing in Flavor Observables 

The F functions that arise from calculating the box diagrams for AF = 2 processes 
are given by [52]: 

7 \ x 2 logXj n ( XiXj 



^ Mw) = (i- Xi )(i- Xj ) v - i XiXj ) + (x.-^a-x,) 2 1 1 - 2Xj + . 

4 log., (l _ 2x XiXj 



{xj - Xi){l - Xj) 2 V ' 4 

(B.l) 

where the corresponding box diagrams have been calculated in Feynman-t'Hooft gauge. 
Since we computed the mass eigenvalues for the top sector at all orders in the v/f 
expansion, we also have the precise values for the F functions. 
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The loop functions appearing in the b — > 37 amplitude are (see, for example, 
Ref. [53]): 



Al(x) 



*o(*) 



3x 3 + 2x 2 -22x 3 + 153x 2 - 159x + 46 
log x H 



2(x- l) 4 
3a; 2 



logx + 



36(x- l) 3 
-5x 3 + 9x 2 - 30x + 8 



2(s-l) 4 & 12(s-l) 3 
The _B S /i + yU~ branching ratio is given by [26] 



Br(5 s //+//-) 



SM 



1 + 



y 



SM 



where 



y 



SM 



y = s i 



ar< — 4 3x t 



+ 



log x t 



x t -l ' (x t - l^ 2 
2 + 2x t - 2x 2 t x t (2 -x t + 2x\) 



log x t 



3 + 2x 



(B-2) 



(B-3) 



8(-l + ar t ) 8(-l + x t ) 2 

Note that these expressions are only valid to order (v/f) 2 . 



logxr + — tan a 
8 8 



.(B.4) 
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